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Abstract 

The random-cluster model, a correlated bond percolation model, unifies a range of important models 
of statistical mechanics in one description, including independent bond percolation, the Potts model 
and uniform spanning trees. By introducing a classification of edges based on their relevance to the 
connectivity we study the stability of clusters in this model. We prove several exact relations for general 
graphs that allow us to derive unambiguously the finite-size scaling behavior of the density of bridges 
and non-bridges. For percolation, we are also able to characterize the point for which clusters become 
maximally fragile and show that it is connected to the concept of the bridge load. Combining our 
exact treatment with further results from conformal field theory, we uncover a surprising behavior of the 
(normalized) variance of the number of (non-)bridges, showing that it diverges in two dimensions below 
the value 4cos^ (tt/VS) = 0.2315891 • • • of the cluster coupling q. Finally, we show that a partial or 
complete pruning of bridges from clusters enables estimates of the backbone fractal dimension that are 
much less encumbered by finite-size corrections than more conventional approaches. 
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1. Introduction 

Percolation is probably the most widely discussed and arguably the simplest model of critical phe¬ 
nomena. Due to a combination of conceptual simplicity and wide applicability which is a signature of a 
problem of very general interest, its many incarnations including bond and site percolation on a lattice 
as well as continuum and non-equilibrium, directed variants, have been the subject of thousands of stud¬ 
ies [1]. In statistical physics and mathematics alike, the quest to understand aspects of the percolation 
problem has led to developments of powerful and beautiful new techniques mm- While the problem 
is well dehned and interesting for any graph or lattice, both finite and infinite, the understanding of 
non-trivial cases is most advanced in two dimensions. There, the scaling limit of critical percolation can 
be related to the Coulomb gas |4] and conformal field theory [5] , leading to exact results for most critical 
exponents and certain correlation functions. More recently, a rigorous approach to conformal field theory 
was pioneered by Schramm who used a mapping introduced by Lowner to construct a way of generating 
conformally invariant fractal random curves, the Stochastic (or Schramm) Lowner Evolutions (SLEs) [5], 
for which a range of properties, including fractal dimensions, can be calculated exactly. In this context. 
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Smirnov and co-workers used the concept of discrete analyticity to establish rigorously that the scaling 
limits of critical percolation [71 [HI and the Ising model |H] on the triangular lattice are indeed conformally 
invariant, and cluster boundaries in these models converge to certain classes of SLE traces. 


The random-cluster (RC) model was suggested by Fortuin and Kasteleyn as a natural extension of 
the (bond) percolation problem, noting that there was a class of models fulfilling the series and parallel 
laws of electrical circuits that also included the Ising model m- Given a graph G = {V,E), it assigns to 
a spanning subgraph {V, A C E) a probability mass (in the following also referred to as RC measure) [TT] 


'^p,q,G[A] 


p\A\{l-p)\E\-\A\^K(A) 

ZKc{p,q,G) 


A e no, 


( 1 . 1 ) 


where K{A) is the number of components and |A| the number of edges in A. The quantity Zrc(Pi G) 
is the partition function of the RC model, corresponding to the sum of unnormalized weights, and VLq 
constitutes the set of all spanning subgraphs or configurations, i.e., Q.a = {A : A C E}. Edges that are in 
A are called open and those in E\A closed. The presence of the cluster-weight factor distinguishes 
(1.1) from the percolation problem and prevents Pp_q_G[^] from being a Bernoulli product measure; only 
for g —>■ 1, where the model reduces to the percolation problem and hence edges become independent, this 
property is restored. Although the cluster weight q can be any non-negative real number, integer values 
of q are particular in that for them the partition function is very closely related to that of the g-state 
Potts model [TH] (see Eq. (2.3) below). For lattice graphs in at least two dimensions, the model undergoes 
a percolation phase transition at a critical value pdq) of the bond probability where, for sufficiently large 
g the transition becomes discontinuous HU. On the square lattice, self-duality allows to deduce the exact 
transition point Pc{q) = Psdiq) = y/q/{^ + y/q) [IS], and the location of the tricritical point is known to 
be gc = 4, beyond which the transition becomes of first order [14] . 


The crucial importance of the RC description for the understanding of critical phenomena is through 
its expression in purely geometrical terms. Hence understanding the geometric structure of the (corre¬ 
lated) percolation problem (1.1) provides a geometric route to the understanding of the thermal phase 
transition of the Potts model. Correspondingly, significant effort has been devoted in particular to inves¬ 
tigations of the structure of the incipient percolating cluster. While initially it was assumed that it was 
a network of nodes connected by essentially one-dimensional links, results regarding the conductivity of 
the critical cluster implied that, instead, the structure is better described by the more elaborate ‘links- 
nodes-blobs’ picture [T5|. If one fixes two distant points A and B on the cluster, those bonds that have 
independent, non-intersecting paths to both A and B form the backbone of the cluster. It is essentially 
the part that would carry current in case a voltage was applied between A and B. The remaining cluster 
mass is in dangling ends. The backbone itself consists of singly-connected or red bonds that destroy 
percolation or, equivalently, destroy conductivity between A and B if cut (the links which meet in the 
nodes) as well as bonds in cycles that are multiply connected (the blobs). These subsets of bonds have 
fractal scaling with associated exponents dp for the cluster mass itself, dsB for the backbone, and dpB 
for the red bonds. Obviously one has dne < dsB < dp, and it turns out that strict inequality holds 
in the generic case and hence asymptotically almost all of the bonds of the incipient percolating cluster 
are in dangling ends, and almost all of the mass of the backbone is in the blobs. In two dimensions, 
exact expressions for dp and djiB are available from Coulomb gas arguments m , but c^bb only known 
numerically [niHH]- 


As we discuss here, the distinction of singly-connected and multiply-connected bonds first suggested 
for understanding the backbone structure [TS] is a more generally useful classification of bonds in a 
cluster configuration. Removing singly connected bonds, which we call bridges^ generates an additional 


^Note that in some of the previous literature on percolation, the term bridge is used to denote the (candidate) red bonds 
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connected component, while this is not the case for multiply connected bonds or non-bridges. Indeed, 
this separation allows us to control the effect of a local edge manipulation (open e-)- closed) on the weight 
of the configuration as in Eq. (1.1|. We will exploit this below to derive expressions for the expected 
number of bridges and related classes of bonds as well as the scale of fluctuations in the former. A further 
partitioning of bridges was recently suggested for percolation in Ref. [18] : branches are bridges for which 
at least one of the connected clusters is a tree, whereas the remaining bridges are dubbed junctions. This 
distinction is related to but not congruent with the links-nodes-blobs picture: most of the red bonds will 
be junctions, but so are some bonds in dangling ends as well. Cutting the branches, on the other hand, 
will remove some part of the dangling ends (the treelike ones), and branches are not typically found in 
the backbon^ As a result of this incongruence, the number of junctions on the percolating cluster grows 
as and not proportional to , where L is the linear size of the system. In contrast to the class of 
red bonds, therefore, both branches and junctions are extensive subsets of the bridges. For uncorrelated 
percolation Xu et al. [m find numerically, however, that removing all bridges, i.e., both branches and 
junctions, leads to tightly connected remainders, corresponding to the blobs, whose size grows as ^ 
As we will show below in Sec. it suffices to remove the junctions only to see the same behavior more 
generally for the RC model. 


As, by definition, the removal of a bridge leads to the generation of a new component and hence 
the breakup of an existing cluster, understanding the properties of bridge bonds is crucial also to the 
understanding of fragmentation phenomena in the framework of a lattice model. This connection was 
investigated in our recent letter [20j . where we studied the fragmentation rate and kernel, and related 
the associated scaling exponents to the more standard critical exponents and fractal dimensions. In the 
present paper, however, we follow a different line of study. In particular, we derive a linear relation 
between the expected density of bridges and the density of bonds; we provide the asymptotic densities 
of bridges in the thermodynamic limit for the square lattice and study the scaling corrections in finite 
systems in general. We employ Monte Carlo simulations based on the Swendsen-Wang-Chayes-Machta 
algorithm CD 123 and a recent new implementation |2il] of Sweeny’s algorithm |24j to confirm these 
results for the square lattice and study the behavior in three dimensions. Building on these results for 
the densities, i.e., the first moment of the bridge distribution, we move on to studying the fluctuations, 
corresponding to the second moment. We reveal an unexpected finite-size behavior leading to a non¬ 
specific-heat variance singularity for q < 4 cos (tt/i/S) , a direct consequence of an intriguing interplay 
between bridges and non-bridges. 

The rest of this paper is organized as follows. In Section we introduce the edge classification, 
derive the bridge-edge identity in general, consider some of its consequences and present a numerical 
study for the special case of the square lattice, together with a detailed analysis of finite-size corrections. 
In Section we use symmetries of the RC model in order to derive connections to other quantities of 
interest in the percolation literature, such as the short range connectivity. In Section]^ we introduce the 
concept of bridge load which allows us to analytically characterize the point of maximal bridge density for 
independent bond percolation. In Section]^ we derive a second-order bridge-edge identity that allows us 
to study the variance of the number of bridges. In Section]^ we compare various new and old methods to 
extract the backbone fractal dimension, with an emphasis on the involved finite-size corrections. Section]^ 
contains our conclusions. 


[161119) . While these are actually backbone bridges, we use the term in a more general sense here to denote singly-connected 
bonds anywhere in the cluster configuration. 

^This is a plausible assumption as the defining property of a branch requires that at least one end of the considered edge 
is a tree, which in the presence of blobs on the backbone is unlikely in the sense of an asymptotically vanishing density. 
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2. Bridge density 


2.1. Edwards-Sokal coupling 

Although in this paper we focus on the RC model itself, its relation to the Potts model is of relevance 
for the physical interpretation of the (geometric) percolation transition in terms of a thermal phase 
transition. In particular, we make reference to how certain observables of the Potts model translate into 
quantities in the RC language. The Potts model describes interacting spins, and assigns one of the q spin 
values in {1,... ,< 7 } to each vertex, attributing a configurational probability or Boltzmann weight 


q,G[o-\ = 


1 


ZpiP,q,G) 


^exp /3 ^ 
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to a spin configuration cr = {ai, ..., cr|y|} G Q = {1, • • • , Here, we write for the indicator 

function, which equals 1 if era, = ay and 0 otherwise. Analogous to Zrc(p, Q, G), the Potts partition 
function Zp(^, q, G) corresponds to the sum of the unnormalized weights. Importantly, in the second 


equality, we need to identify p = 1 — e The relation between Eqs. (1.1) and (2.1) is established by 


augmenting phase space to include both, spin variables cr and subgraph variables A, leading to the joint 
probability 


^ P.9 


ck. A] = 
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that is also known as Edwards-Sokal coupling [25] . It has the crucial property that the marginal distribu¬ 


tion on the spins alone reduces to the Potts weight (2.1), while the marginal on the subgraphs gives the RC 


weight Similarly, one can derive the conditional probabilities for the spins given the bonds or vice 

versa m- As was already shown by Fortuin and Kasteleyn [10] . the corresponding partition functions 
coincide up to a (trivial) multiplicative factor if we relate the parameters /? and p via p = 1 — e~^, 


Znc{l-e-^q,G) =e-f^\^\ZpiP,q,G) for g G {2,3, • • • }• 


(2.3) 


This coupling is the basis for the Swendsen-Wang cluster algorithm [21] and its generalization to non¬ 
integer q due to Chayes and Machta [22]. The above coupling and the identity (2.3) allow in particular 


to relate expectation values in the two ensembles. For instance, the internal energy in the Potts model is 
up to constants identical to the expected number of edges in the RC language. 


Up,q,G = 


1 d\nZp{f3,q,G) 

'M ^ 


-^|^Ep,,,g[A/'], for p = 1 - 6 ” 


(2.4) 


where A/’(A) = |A|/|A| is the proportion of open edges and Ep^q^c]-] denotes the expectation with respect 


to the RC measure (1.1). Similarly, higher moments of the energy distribution such as the specific heat as 


well as magnetic observables of the Potts model can be expressed in terms of expectation values in the RC 
language m- Due to the conditional measures, correlation functions can also be related. Consider, in 
particular, the event x ga y, that is the existence of a path in A connecting x and y. From the conditional 


measures derived from Eq. (2.2) it is easy to see that 


Ep,q,G[l - t{a^=ay}\A] = 


9-1 


1 - 1 


jxAy} 


(2.5) 


That is two vertices can only have unequal spin when they are not connected in the bond configuration 
and were assigned two different (random) colors. Summing over the edges shows that the nearest-neighbor 
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connectivity defined as the sum 


£(A) = — V Ir A 1 
^ \E\ \x^y] 

' ' {x,y)&E ^ J 


( 2 . 6 ) 


is essentially the energy [IIlEelET] 

\E\ q-l(\E\ 


Ul3,q,G = 


fOTp=l-, 


q 

Clearly for hypercubic lattices |ill|/|y| is bounded (more precisely equals d, where d is the dimensionality). 


2.2. Edge classification 


Intuitively, (1.11 indicates that for g > 1 (g < 1) the model favors configurations with a larger (smaller) 
number of components, compared to the case g = 1. In this paper we investigate how the cluster structure 
is influenced by the cluster weight g, in particular with respect to the fragility of clusters. The relevance 
of individual edges for the connectivity of A C E depends on whether they are pivotal or non-pivotal. 
According to our definition, pivotal edges change the number of components upon removal/insertion, i.e., 
for a pivotal edge e we have K(A^) = K{Ae) — 1, where A® (Ag) is the configuration obtained from A by 
opening (closing) e (note that necessarily one of A® and Ae is equal to A). Hence for a pivotal edge there 
is no alternative path connecting the vertices incident to it. As discussed above, for the purposes of this 
paper open pivotal edges are denoted as bridges, while we call pivotal closed edges candidate-bridges. 
Likewise, we denote open (closed) non-pivotal edges as non-bridges (candidate-non-bridges). Let B, 
C, B, C denote the set of bridges, non-bridges, candidate-bridges and candidate-non-bridges, respectively. 
We stress that these sets depend explicitly on the configuration A. We denote the corresponding densities 
as H = |i?|/|A|, C = \C\/\E\, B = \B\/\El and C = |C|/|A|. 


2.3. Derivation of the bridge-edge formula 


In this sub-section we show how, for arbitrary (finite) graphs, the expected densities of all the above 
mentioned edge-types B, C, B, and C can be related to the expectation of the edge density. We establish 
these identities by utilizing a differential equation known as Russo-Margulis formula [2H1[2S], which applies 
to expectations with respect to product measures such as the probability density of the bond percolation 
model, corresponding to the g = 1 RC model. The formula intuitively quantifies the “response” to an 
infinitesimal change in the bond density. Using this approach to quantify the response in the number of 
connected components K allows us to study the bridge density. We then proceed and apply the idea to 
the RC model with generic g > 0, by exploiting the fact that the partition function of the model can be 
interpreted as a particular expectation in the percolation model. 


First consider the percolation model, i.e., the measure (1.1) with g 
tions as In this setting, the Russo-Margulis formula reads 


1 , for which we denote expecta- 


^Ep,G[A] = ^Ep,G[5eA], 
^ e^E 


(2.7) 


where X{A) : H — >■ K is an arbitrary observable and the quantity SeX is called the influence of e on A 
(or the discrete e-derivative of A) and is given by 


((5eA)(A) = A(A^)-A(Ae). 
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The graph specific geometric properties of the model under study are “encoded” in the influences. In 
order to simplify the subsequent analysis, we show in Appendix the following bijection identity: 

Ep^alSeX] = -Yl [^(^) - ■ ( 2 . 8 ) 

P ACE; 
eC A 


We now apply this identity to the cluster-number observable K, for which K{A) — K{A^) = —1 if e is 
a bridge, e S B{A), and 0 otherwise. Application of the Russo-Margulis formula (2.7) and Eq. (2.8) to 
X = K then yields 


^Ep.G[if] = 5] Ep^G[5eK] = - 51 E [K{A) - K{A,)\ = - - E (^B\ = 

P eeE P eeE aqe, P geE P 

ee A 

(2.9) 

This differential equation allows us to determine the cluster numbers Ep ^[Ar] once the p-dependence of 
Ep_g[;B] is known. We remark that this relation was derived before in [3D], however with a different target 
application in mind. As already mentioned, we can extend the above to the RC model and obtain more 
generally the bridge density Ep^^j^cl'B] of the RC model itself. To start with, we note that the RC model 
partition function is ZYic{p,q,G) = Ep^cl?^]- In other words, the RC model partition function is the 
expected value of in the bond percolation model with parameter p. Hence we can apply the Russo- 
Margulis formalism to q^. Using the same reasoning as that for Ep^G[Ar] above, we show in Appendix 
that 

O -I 

— \ogZnc{p,q,G) = |a;|^Ep.,.g[H]. (2.10) 

On the other hand, direct differentiation of the partition function Zrc(Pj G) yields 


_5 

dp 


logZRc(p,g,G) = Qep,,,g[A /1 - . 


( 2 . 11 ) 


Comparing both expressions ( |2.10| ) and (2.11), we derive the following relationship for the RC model 
(referred to as bridge-edge formula). 


E. 


'P,9,G 


[B] = , 


E. 


•p,q,G 


[N]-p 


{l-p){l-qy 


( 2 . 12 ) 


which is of central importance for the following analysis (see m for some related identities developed 
along different lines). We remark again that this holds for any graph and for any p G (0,1) and q G (0, oo)\ 
1 (the cases 9 —>■ 0, 1 are discussed below). Due to the obvious identity Ep_q_G[Af] = Ep^q^G[B\ -I-Ep^^^Gp], 
one immediately arrives at an analogous equation for the density of non-bridges. 


E, 


P,q,G 


[C] = 


ipq-p- g)Ep,,;,G[A /1 -p 
(l-p)(l-q) 


(2.13) 


We remark that in principle one can use the Russo-Margulis formalism to study the (partial) p-derivative 
of expectations of other observables in the RC model: 


—Ep,,,g[a:] 


^ / Ep,G[g^X] \ 

dp V Ep,G[g^] ) 


^p,g[q^] 


9-1 


Ep,g,G[-A]Ep^,^G[|S|], 


which shows the relevance of bridges in general for the RC model, 
direct consequences and applications of ( 2 . 12 ): 


Furthermore, let us consider a few 
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1. Firstly, notice that (2.12) recovers the correct results for the limit of uniform spanning trees, where 
all open edges are bridges. This tree model can be obtained in the RC model in the limit q,p ^ 0 
for q/p —>■ 0, which for the square lattice includes the critical line Psdig) [II]- 

2. It is possible to establish an upper bound for ¥.p^q^GW\ for g > 1 and any p on any graph. Using 
general comparison inequalities for the RC measure |32j . it is possible to show that 


Ep,( 3 .g[A /1 > ]Ep(p,q),l,G[A/'] = 


{l-p)q+p 


= Pip,q), 


where p{p, q) < p ior q > 1. Together with (2.12) this yields the desired upper bound for ¥.p^q^Q[B] 

p - p{p, q) 


E. 




(l-p)( 9 -l) {l-p){q-l) 


= pip,q)- 


As is shown below for the Ising case q = 2 with an exact solution this bound is tight in the limit 
p —>■ 0, but clearly not so for p > psd- 


3. Another interesting consequence follows when one recasts (2.12) to express the edge density in terms 
of the bridge density 

^p,q,GW] =P - (l- 1)(1 -p)Ep.g,G[^]- (2.14) 

This emphasizes the importance of bridges for q ^ 1 as the source of finite-size corrections in 
]Ep,(j,g[AC], a manifestation of the introduction of correlations (deviation from the product measure). 
Moreover it also explains, in the setting of the RC model, how thermal fluctuations in the number 
of edges vanish in the percolation limit q —>■ 1, namely by a vanishing factor (g — 1). This has been a 


plausible assumption in the literature so far, but to our knowledge, relation (2.14) is the first exact 
statement, see for instance 


2.4- The square lattice 

For the infinite square lattice the self-dual point Psd{q) = ^/q/{^ + \/q) is criticaQ For the Potts 
model at criticality, a mapping to an exactly-solved six-vertex model gives extensive analytical results 
[Hdl]; it is assumed that by analytic continuation these carry over to general values of q. In particular. 


the critical internal energy density in this case is = 1 -|- ^jyfq and hence from Eq. (2.4) one 

deduces that = 1/2, where we defined in general E^^gI’] = ^Pcfg),?!’] graph G (analogous 

for probabilities), where of course the critical line is given by Pc{q) = Psd{q) for Z^. It follows from (2.12) 


and (2.13) that one has 


Eg,Z2[;B] = 


2 ( 1 + V9)’ 


Eg,Z2 [C] = 




2 (1 + V9)' 


We note that Eqs. (2.15) can be written in terms oi Psdiq) and Eg^z 2 [A/], as 


Eg,z2 [13] = [1 - Psd (<?)] Eg_z 2 [A/"], Eg_z2 [C] = Psd {q)^q,z^ [A/]. 


(2.15) 


(2.16) 


The above expressions allow us interpret Psd{q) as the expected fraction of non-bridges among all open 
edges. Clearly, the remaining fraction 1 — Psd{q) of open edges have to be bridges. We note that 
because ofeeR=>eSAwe have Pp,g,G[e & B\e G A] = Pp,g,G[e G il]/Pp_g_G[e G A] and hence in 
particular Pg^z 2 [e G B\e G A] = 1 — psd{q) and similarly Pg^z^ie G Cje G A] = Psd{q)- Interestingly, an 


®Note however that albeit being numerically and physically well supported, this has only recently been rigorously 
confirmed for the case q > 1 m- 
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immediate consequence of (2.151 is that our RC result reduces in the percolation limit, g —>■ 1, to ^2 [B], 
—>-1/4, which recovers a result recently derived in Ref. [T5] . 

Going beyond criticality, for the Ising model, corresponding to g = 2, we can use the available exact 


solution for the internal energy and the relation (2.4) to find the exact p-dependence of the densities 


of edges as well as bridges and non-bridges. This allows us to study, in an exact setting, the fragility 
as the bond-density/temperature is varied. As the corresponding expression for the internal energy is 
not very instructive we refrain from reproducing it here, see, e.g.. Ref. [14] for details. Instead, we 
show in the upper panel of Figure the exact asymptotic results for iF'. For comparison, we also show 
numerically estimated values of the density of bridges for other cluster weights (for details regarding the 
numerical procedure, see Sec. 2.6 below) in the lower panel of Fig. As it is clear from Fig. the 
density of bridges is not maximized at the critical point but for some p{{q) < Psd)?)- We shall give an 
explanation for this phenomenon in terms of the relationship between nearest-neighbor connectivity and 
bridge density for the percolation case g = 1 in Sec. |2.5[ Regarding the results for the Ising model, we 
note that it is also possible to produce exact expressions for the internal energy on finite lattices [M] , 
resulting in corresponding exact results for the bridge densities on Z|, the L x L square lattice with 
periodic boundary conditions. 

It is apparent from the lower panel of Fig. that a characteristic of the critical point is the extremal 
slope of Ep q xi [;B]. Moreover, the relation \2.12\ suggests that the slope becomes singular due to a 
thermal singularity for q>2. Indeed, the occurrence of a singular slope can be explained by using (2.10) 
connecting Mp^q^c[B] to the first p-derivative of the log (Zp^q), which allows one to derive 


dpEp^q^a[B] = 


p5plog[ZRc(p,g,G)] 1 


\E\{l-q) 


-Ep,g,G[^]- 

P 


(2.17) 




pa(q) 


this latter relation together with (2.4) yields the asymptotic scaling 9pEp ^ ['®]lp=Pad( 9 ) 

logL for q = 2 [T3] and for all other values in g G (0,4], with multiplicative logarithmic correction 
proportional to log“^/^ (L) for the tricritical case q = 4 PSII55H57] . Here, / is the free-energy density. We 
note that Coulomb gas arguments yield expressions for ajv for q € (0,4] that predict a/v >0 whenever 
g > 2 |4], implying a finite slope for g < 2 even in the limit L —>■ oo. We refer the reader to Section |2.7 
for more details about size dependent effects and the universality of the above arguments. 


2.5. The percolation case 


Returning to the case of general graphs, we see that the bridge-density formula (2.12) is singular for 
g = 1. To deduce the correct result in the percolation limit g —)■ 1, we use L’Hopital’s rule to find 


limEp^q^clB] = 


1 

1 — p g-^-i ~ 1 


-T— 

l-pq^ioq 


,.G[An = -j^Covp^g[K,N]. (2.18) 


-P 


•p,q,G[Ff] —?► p as g —>■ 1. The l.h.s. being a density, which in particular must be 


Here, we have used that E 
non-negative, shows that the covariance between the number of components K and the density of edges 
Af is non-positive. This is plausible, because adding edges can never increase the number of components. 
In fact on more mathematical grounds, for the monotone case, g > 1, this can even be proven, as it 
follows from the Fortuin-Kasteleyn-Ginibre (FKG) inequality [TT], applicable to the RC model whenever 
g > 1. Applied to K and Af, the FKG inequality yields Ep^q^G[Kff] < Ep_q_G[A']Ep_q_G[A/], which in turn 
particularly applied to g = 1 shows that CoVp^GiKAf] = Ep^G[KAf\ — Ep^G[Ar]Ep^G[A/] < 0. Note that the 
covariance in (2.18) was analyzed for the percolation model in Ref. |38] . however without establishing a 


connection to the bridge-density. There, the authors discuss the size dependence of CoVp^^^d [K,Af], where 


















Figure 1: Upper panel: Exact densities of edges, bridges and non-bridges for the q = 2 RC model on The red line shows 
the bridge density, the blue line denotes the non-bridge density and the green curve represents the edge density. The inset 
shows the expected fraction of bridges (red) and non-bridges (blue) among all open edges. Lower panel: For comparison 
we also show numerical estimates for the bridge density on Z|^ for cluster weights in [0,4] and L = 64 whenever q ^ 1 and 
L = 2048 for q = 1. The vertical lines show the location of the self-dual (critical) points for the relevant values of q. 


zi is the d dimensional hypercubic lattice with linear dimension L and periodic boundary conditions. 


we provide an explanation for this, in terms of the bridge density. 

Before we proceed to discussing the behavior on finite lattices in the RC setting, we derive a symmetry 
relation for the density of bridges, valid for Z^, which for the special case p = 1/2 implies the previously 
shown result that Ei/2,z2['S] = 1/4 for critical percolation on Z^. Let us naturally view the vertex set of 
Z^, denoted by V (Z^), as the set of (ordered) tuples, whose entries are integers, and the set of edges of Z^, 
denoted by E(Z‘^), is the set of pairs of vertices which differ in precisely one coordinate by ±1. Then we 
define a sub-box of Z^, denoted by G„ = (14, i?„) for n > 1, where 14 = {{x, y) G V(Z^)| —n < x,y <n} 
and En = E{Z‘^) n l4. This allows us to consider a sequence of sub-boxes (in n), that “approaches” 
Z^. It is not hard to see that each of the G„’s is planar, hence we can use Euler’s formula for planar 
graphs, stating that for any A C En, K{A) = |14| — |4| -|- F{A) — 1, where E{A) is the number of faces 


The authors find that the quantity has a leading finite-size correction that scales as In Sec. 2.7 
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Figure 2: A graph (left panel) with one of its depth-first search trees (right panel). Red (hollow) and blue (solid) edges 
correspond to bridges and non-bridges, respectively. The dashed edges in the right panel are back-edges. Note that the 
depth-first search tree together with its back-edges can be decomposed into two chains, corresponding to the two back-edges, 
and any non-bridge is in at least one chain. 


in any planar embedding of the graph {Vn,A). Taking expectations on both sites and dividing by \En\, 
we obtain 

Ep^ajK] _ |K| 1 

\En\ \E^\ ^ \En\ |S„r 

Furthermore it can be verified that for any planar graph the number of components in the primal graph 
{V^A) equals the number of faces induced by the dual configuration A* in the dual graph {V*,A*), in 
other words we have applied to this case, Fq^{A) = iFc* (A*), see e.g. [32] • Thus, we can re-write the 
above 

Ep^cAK] _ |K| 1 

\En\ \En\ ^ \EA \En\' 


which uses the duality of bond percolation, i.e., Pp_G(^) = Ei-p,G*(^*)i valid for any planar graph G. 
We can now differentiate both sides with respect to p and apply the Russo-Margulis formula (2.91 to 
obtain 

=P - ^ ^^ Ei_p^g* [g]. 


Now assuming that lim„_>oo Ep^c^ [B] = Ep ^2 \B] and exploiting self-duality of 1? in the form of the 
plausible assumption that lim„_>.oo Ei_p,G* [^8] = Ei_p_z 2 [;B], we obtain 


Ep,z2[^] =p - ^ ^^ Ei_p^z2[g]. 

Finally, for the particular choice of p = 1/2 we recover the (asymptotic) result Ei/ 2 ^ 7 ? [B] = 1/4 HHj. 
Moreover the above duality result implies a symmetry of Ep ^2 [B] around the self-dual or critical point 
p = 1/2. In other words determining Ep Z 2 [S] on [0,1/2] suffices to obtain Ep Z 2 [,B] on the entire interval 
[ 0 , 1 ]. 


2. 6 . Numerical analysis for finite lattices 


In order to confirm the asymptotic results (2.15) for the square lattice and to study the finite-size 
corrections for finite lattices, we performed Monte Carlo simulations of the RC model on the Lx L square 
lattice with periodic boundary conditions, in the range of continuous phase transitions 0 < g < 4. We 
used a recent implementation [23] of (the Metropolis variant of) Sweeny’s algorithm [23] for g < 1 and 
the Chayes-Machta-Swendsen-Wang algorithm [22] for q> 1. We determined the number of bridges in 
a given configuration by means of the algorithm introduced in Ref. [39]. In contrast to the approach 
used in [nmn], this is a linear time algorithm applicable to any finite graph. It does not depend on a 
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Figure 3: Critical bridge- and non-bridge density for a L = 32 system shown as square symbols. The circles show the 
asymptotic values derived from a finite-size scaling analysis using the ansatz ( |2.20[ |. For comparison the solid (blue) and 
dotted (red) line show the exact asymptotic values of non-bridges and bridges, respectively, according to l |2.15| |. 


medial graph analysis, which facilitates the study of higher dimensional and non-planar systems. Here 
we briefly describe the algorithm and refer the reader to Ref. [39] for more details. The main idea is to 
construct a depth-first search (DFS) tree for any component in the spanning subgraph {V,A). Any edge 
in A that is not in a DFS tree is called a back-edge (there are precisely |A| — |H| + K{A) back-edges). 
Because a back-edge connects vertices already connected in a DFS tree, it follows that any back-edge is 
a non-bridge in (H, A); see Fig. [^for an illustration of this situation. In order to be able to determine 
whether a given edge in a DFS tree is a bridge in (V, A), the algorithm introduces a chain decomposition 
of the graph, where for any back-edge there is exactly one chain (ignoring chains consisting only of one 
vertex). Now, if a given tree edge is part of such a chain, then it is in at least one cycle and hence a 
non-bridge. Finally, a careful construction that avoids the iteration over overlapping chains ensures that 
the algorithm classifies all edges into bridges and non-bridges in linear time. 


2.1. Finite-size corrections 


The bridge-edge relation (2.121 also allows us to understand the finite-size corrections to the bridge 
density. It relates to the density of edges which, in turn, is related to the internal energy 


via (2.4). Finite-size corrections to the energy density U/ 3 ,q,G close to a point of a second order phase 
transition, on the other hand, are widely studied and well understood from the theory of finite-size scaling. 
To formulate it, we need to associate a linear size L to G, such that |V| = as is the case for a lattice 
graph with spatial dimension d. Then, the standard ansatz for the singular part fs of the free-energy 
density / = — log [Zp(/3, g, G)]/(/3L‘^) as a function of the reduced temperature t = {T — Tc)/Tc, the 
external field h, and the leading irrelevant field v is given by [41] 


fs{t, h, L) = L-'^F{tL^/'', . 


(2.19) 


The internal energy ^ gd can be expressed in terms of the first derivative of / with respect to t, 
evaluated at t = h = 0. Furthermore it is plausible to assume mma that the non-singular part of / has 
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no size-dependence and directly yields the value obtained in the infinite-volume limit. We hence expect 
the following scaling of the bridge density: 




( 2 . 20 ) 


Here, we introduced the L-subscript to denote the expectation for the model on any lattice graph with 
linear dimension L, corresponding to the universality class of the g-state Potts model on Z^. For the 
case of Z|, the asymptotic density E,j_oo[^] is the one of Z^, as given by Eq. (2.15). For other 2D 
lattices such as the triangular and honeycomb lattice Eq^co[B] can be derived in a similar way using 
well-known expressions for the critical internal energy there, see, e.g., Ref. | 
itself, however, is universally valid. 


The scaling form (2.20) 


We now turn to a more detailed analysis of the corrections for the case of the square lattice. We used 
the above bridge detection algorithm to numerically determine the size dependence of the bridge density 
for the critical RC model on Z| for a number of values of the cluster weight q. Testing for scaling of the 
form (2.20), we fitted the ansatz a + bL'^ to our Monte-Carlo estimates, using the method of least squares. 
Up to cluster weights of about q « 3.25, all results were found to be in agreement with the implication 
of Eqs. (2.20) and (2.15) that a = 1/[2(1 -I- ^/q)] and c = Xjv — d, where d = 2. Fixing the parameter 
a at this exact value, yields more precise fit results for the exponent c. These data are shown in Fig. 
together with the exact values oi 1/v — d known from the Coulomb gas mapping [U [16]. The deviations 
observed for q > 3.25 are attributed to the presence of higher-order corrections to the form (2.20) which 
increase in strength on approaching the tricritical point q^ = 4. 


Regarding the prefactor b of the leading term extracted from the fits, we find that it is 

negative for the bridge density, hence the asymptotic result 1/[2(1 -I- ^)] is approached from below. 
The corresponding constant for [B] is positive such that the limit is approached from above. This 
holds for all values analyzed in q G [0,4]. Thus in illustrative words, the configurations on Z| still 
“feel” the “extra” edges closing the lattice to form a torus and hence they are typically less “fragile” 
than corresponding configurations on 1?. Interestingly, when one considers the finite-size dependence of 
[-^l one finds an excess of edges (relative to the asymptotic result for Z^) for g > 1 and a shortfall 
for g < 1. Now, due to the fact that Ep_q_( 3 [A/] = Ep^q^Q[B\ + Ep^q^Q[C], we have that in particular the 
amplitudes of the finite-size corrections of E^ ^2 \B] and E^ [B] cancel each other for g = 1, whereas the 
modulus of the amplitude for Eq j 2 ^[B] is larger (smaller) than the one for Eg_ 2 |[C] for g < 1 (g > 1), 
respectively. We omit the details on the size dependence of E^ 2 |[C]i ns no new mechanism appears and 
the above results on E^ [B\ can be easily adapted to this case. 

Following Xu et al. CHI, we extend our study of finite-size corrections on the torus Z|^ from bridges 
to the larger class of “type-1 edges”, that is edges that that have both of their associated loop arcs in the 
same loop in the associated loop configuration. While on a planar lattice all such bonds are bridges, this 
is not the case for Z|^, where also non-bridges can be of type 1. Such edges are hence called pseudo¬ 
bridges [IHl SO]. If we denote by Li(A) and L2{A) the set of (open) edges in (14,^4) that have both 
of their associated loop arcs in the same loop (type-1 edges) or in two different loops (type-2 edges), 
respectively, we have that 


^P,q,zl “ ^p,q,zl [e G Li,e G B] + ]Pp,5,z| [e G Li,e G C] — [e G B] + IPp,9,z| [e G Li,e G C]. 


( 2 . 21 ) 

The second equality follows from B C Li, that is any bridge has necessarily both loop arcs in the same 
loop (as it cannot enclose any face). Equivalently, we have Ep^^g 2 = Ep^q^^ 2 ^[B] + Ep^z 2 [7^], where 
V{A) = \e G E ■. e G Li{A),e G C'( 24 )| /m, i.e., the density of pseudo-bridges in A, and £i denotes the 
fraction of edges that are of type 1. In order to understand the general size dependence of E^ 2 |[^i], 
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Figure 4: Comparison of numerically extracted exponents of the leading finite-size correction for the density of bridges and 
pseudo-bridges. The solid and dotted line show the exact Coulomb-gas value for pseudo-bridges and bridges, respectively. 
The deviations for values of q around 4 are caused by strong sub-leading finite-size corrections. 


determined the number of type -1 edges in our numerical simulations and fitted a finite scaling ansatz 
Eg ^2 [ii\ = a + + cL~'^^ to our Monte-Carlo data. The resulting estimates are consistent with 

a = 1/[2(1 -I- y^)] and ei = 2 — l/i/ as well as £2 = X 2 - In a second step, we also performed fits with c 
fixed to its exact value and ei and £2 fixed to the values implied by the Coulomb gas for the identihcations 
£1 = 2 — 1/v and £2 = X 2 , such that b and c were the only free parameters in a now linear fit. As it is 
apparent from the parameters collected in Table the quality of these fits is very good, even including 
the smallest system sizes L = 4, 8 , 16, 32, strongly suggesting the following asymptotic form for the 
square lattice 

[i,] ~ E,,z 2 [B] + + cL-"=, ( 2 . 22 ) 

where again = 1/[2(1 + \/q)]- Figure[^shows our numerical data for the two cluster weights 0.5 

and 2 , together with the corresponding best hts. 


It is interesting to compare the general form (2.22) to the results found in Ref. [18] for the percolation 
case, corresponding to g —>■ 1. There, finite-size corrections to the bridge density on the torus were 
found to decay with exponent —X 2 = —5/4. While this is numerically consistent with our hndings since 
l/i^ — d = —X 2 for percolation [TTl IUl - 146] . it is clear from the present analysis that the relevant exponent 
describing finite-size corrections to the bridge density is 1/v — d. Indeed, the values of —X 2 and 1/v — d 
strongly differ for q ^ 1, cf. Fig. Moreover, for the square lattice with periodic boundary conditions 
Xu et al. [TH] showed rigorously that Pi/ 2 ,i,z|[e € Li] = Ei/ 2 ,i,z| [I^i] = 1/4, independent of L. This is 
consistent with our general form ( 2 . 22 ) for the percolation case with 1 /v — d = —X 2 if the amplitudes 
b and c cancel. Indeed, this is what we find from the fit data in Table which clearly support the 
statement that 6 = —c as <; —>■ 1, as required by the rigorous finding that "^\j 2 .\,T? [^ 1 ] = 1/4. Hence, it is 
a particular cancellation of finite-size effects that leads to the size-independence observed for percolation. 
Further, it is known that for g —>■ 0 with p = psdiq), one recovers the uniform spanning tree model, for 
which clearly no pseudo-bridges exist, that is one can easily verify that, on £^^^2 [^l\ —t 5 — \L~'^ 
for g —0, which is in agreement with (2.22), due to 1/^ —>■ 0 and the numerical observation (see Table 
that a —>■ —1/2 and 5 —>• 0 for g —>• 0. We remark that the htting estimates for g close to one need to 
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q 

b 

c 

X^/Nd.o.{. 

-^min 

Q 

3.5 

-0.1410(7) 

0.504(7) 

1.0354 

4 

0.4135 

3.5 

-0.1405(8) 

0.49(2) 

0.9722 

8 

0.4730 

3.5 

-0.1406(9) 

0.50(4) 

1.0155 

16 

0.4270 

3.5 

-0.140(1) 

0.4(1) 

1.1116 

32 

0.3514 

2 

-0.218(3) 

0.393(7) 

0.5554 

4 

0.9007 

2 

-0.219(3) 

0.40(1) 

0.6303 

8 

0.8183 

2 

-0.219(4) 

0.40(2) 

0.7086 

16 

0.7173 

2 

-0.218(6) 

0.38(4) 

0.8504 

32 

0.5580 

1.25 

-0.27(1) 

0.32(2) 

0.8451 

4 

0.6199 

1.25 

-0.26(2) 

0.31(3) 

0.9229 

8 

0.9229 

1.25 

-0.25(3) 

0.30(4) 

0.8482 

16 

0.5819 

1.25 

-0.31(4) 

0.41(7) 

0.4951 

32 

0.8607 

0.9 

-0.21(3) 

0.19(3) 

0.4983 

4 

0.9356 

0.9 

-0.22(5) 

0.20(4) 

0.5010 

8 

0.9155 

0.9 

-0.23(7) 

0 .21(6) 

0.5237 

16 

0.8748 

0.9 

-0.3(1) 

0.2(9) 

0.6113 

32 

0.7692 

0.5 

-0.336(8) 

0.192(4) 

1.0286 

4 

0.4203 

0.5 

-0.32(1) 

0.188(5) 

1.0466 

8 

0.4019 

0.5 

-0.32(2) 

0.188(7) 

1.2546 

16 

0.2502 

0.5 

-0.33(4) 

0.190(9) 

1.5533 

32 

0.1332 

0.05 

-0.429(6) 

0.033(1) 

0.7662 

4 

0.7073 

0.05 

-0.41(1) 

0.033(1) 

0.6955 

8 

0.7575 

0.05 

-0.43(3) 

0.033(1) 

0.6882 

16 

0.7365 

0.05 

-0.43(7) 

0.033(1) 

0.8571 

32 

0.5522 

0.005 

-0.466(5) 

0.0040(4) 

0.9548 

4 

0.4979 

0.005 

-0.49(1) 

0.0045(5) 

0.6865 

8 

0.7663 

0.005 

-0.49(3) 

0.0046(6) 

0.6857 

16 

0.7389 

0.005 

-0.38(8) 

0.0041(7) 

0.5806 

32 

0.7948 


Table 1: Estimates of the fitting parameters b and c in the finite-size scaling ansatz ( |2.22[ | of the density of type-1 edges. 
Q denotes the quality of fit m or confidence level, that is the probability that would exceed the empirical value, under 
the assumption that the imposed statistical model is correct. In what follows we refer to fits with confidence level > 10% 
simply as “good” fits. 


be treated with caution, as here the two exponents Xjv — 2 and —X 2 come very close to each other, and 
hence it is numerically very difficult to distinguish the two contributions and the corresponding constants 
a, b. However, outside a suitable “safety”-window around g = 1, it appears that both a, b are increasing 
with q. 

We close the analysis of size-dependent effects by establishing a relationship for percolation between 
pseudo-bridges and a particular arm event which shows the observed decay. A pseudo-bridge is a 

non-bridge that resides on a cross cluster that winds simultaneously along both directions on the torus 
[18| . Furthermore, any pseudo-bridge has both associated loop-arcs in the medial graph in the same 
loop. In other words, a pseudo-bridge is a non-bridge on a cross cluster that is pivotal to the existence 
of the cross cluster, i.e., if A is a configuration that contains a cross cluster, then removing a pseudo¬ 
bridge e implies that Ag has no remaining cross cluster. Therefore pseudo-bridges are in a certain sense 
the “bridges” of the cross cluster. We can translate this into a (polychromatic) four-arm event [Sj |47j 
as follows: Any edge, say e, that is pivotal to the cross cluster has the property that starting from e 
two paths of open edges extend to a distance L/2 (here we use the toroidal geometry and in particular 
the absence of boundaries). Since e is pivotal for the cross configuration, there must also be two paths 
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Figure 5: Estimated size dependent deviation ^2 [^ 1 ] — l/[2(y^+ 1)] for q = 0.5 and q = 2 for the self-dual RC model on 
together with the best fit to the scaling form corresponding to (2.22 l All available system sizes starting from L = 4 
were included in the two fits and yielded X^/-^d.o.f. = 1.0286 and 0.5554, for q = 0.5 and 2 respectively, where A^d.o.f. = 14, 
is the number of degrees of freedom for both cluster weights considered. 


of closed (dual) edges starting at the dual vertices associated to e and extending to a distance L/ 2 , 
which precisely ensure that there can be no alternative, e-avoiding, path that would preserve the cross 
configuration property upon removal of e. Thus what we constructed is a particular instance of an poly¬ 
chromatic four-arm event, of alternating primal and dual “color”, in the annulus centered around the 
edge e with outer radius L/2. The probability of such an event decays for large L as L , where x\ 
is the poly-chromatic four-arm exponent, which equals 5/4 for critical percolatiorj^ Using the fact that 
0 : 4 ^^ = X 2 ^ where the latter is the two-arm exponent as defined for instance in |46j and referred to by Xu 
et al. in ng we have the appearance of the exponent X 2 - 

We note that the notion of pseudo-bridges and the corresponding scaling analysis given above is 
specific to the case of locally planar lattices. The concentration on the square lattice was for practical 
reasons only, however, and we expect the same results, only with different amplitudes, for the cases of 
other lattices such as the triangular or honeycomb cases with periodic boundary conditions. 


3. Other types of pivotal edges 

So far our analysis has focused on open edges, that is bridges and non-bridges. Surprisingly, the 
analysis of closed edges turns out to be very interesting too, in that it allows us to link several seemingly 
unrelated quantities, studied in the literature in different contexts, to the density of bridges. Moreover 
the study of closed edges provides a further probe for the cluster structure, for instance a large density 
of candidate non-bridges suggests that the clusters are more likely to self-entangle than to overlap with 
other clusters. In order to analyze the corresponding densities consider a given closed edge e = {x,y). 
Clearly we have that when x y then e is a candidate bridge (and vice versa), i.e., a candidate bridge 


^In fact the value can be rigorously established for critical site percolation on the triangular grid, thanks to the celebrated 
work of Smirnov [8]. 
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corresponds to a pair of nearest-neighbor vertices that are not connected. This observation allows us to 
state the following almost trivial but useful identity that holds for arbitrary e = {x,y) G E: 

Pp,g,G[e G B]= Pp, 9 ,G[a^ ^ y]- 


In fact, it is possible to relate the r.h.s. above to Pp,g,G[e G B], the probability that e is a bridge. To see 
this, note that for fixed e = {x, y) any configuration ACE that has the property x y can be related 
one-to-one to the configuration A + e, in which e is a bridge. Furthermore, as the insertion of e reduces 
the number of components by one, it follows that Pp^q_G[^] = <z(l ~ Py^p,q,G[^‘^\/p which in turn implies 


9(1 -P) 


E' 

ACE 


Pp,(J,g[^ + ^]^{eeB(A+e)} — ^ -^Pp,(3,G[e G B], 

_ (3.1) 


Pp.9.G[a; ^ y] = E r A^I ^ 

ACE 1 1 

Hence we have Pp,q_G[e G B] = ‘^p,q 

the definition (2.6) of the nearest-neighbor connectivity £ we obtain the general result 

"Ep.g.G[Af] 


,G[e G B], and therefore using the bridge formula (2.12) and 


E. 


= 1 - K,q,G{£] = 


The corresponding asymptotic values for the square lattice follow immediately: 

V? „ 1 


[B] = 




Eg,Z2 [C] = 


2(1+ V9)' 


(3.2) 


(3.3) 


As a result of the relation to the edge density, we again conclude that the leading finite-size corrections 
are given by These findings are in line with those of Hu et al. in Ref. m, who find for percolation 

that for L —>■ 00 on one has Ep ^a [£] —>■ 3/4. On the other hand, it was shown in Ref. [IH] that the 
asymptotic bridge density for critical percolation on Z| is 1/4 (which also follows from our result (2.12)). 
These two results are hence clearly consistent with the relation ( |3.1[ ). Finally from (3.2) and (3.3) one 
concludes that 

Eq ['^^1 —>■ . ^ ' for L —>• 00 , 

' 2(1 + ^) 


which is consistent with the results of Ref. m- 


4 . Maximal bridge density for percolation 

It is apparent from Fig. that the expected bridge density is maximal for pf(g) < Psd{q) 

for any q G (0,4] [20] • The reason for this is a priori not clear and here we focus on the percolation case 
<7=1 and obtain an expression for ^Ep gI-B], that is for general graphs G, to determine its zero(s) one 
of which must correspond to the location of the sought maximum. In principle we could again apply the 
Russo-Margulis formula directly to jHj to obtain the required derivative. However, there is an alternative, 
more geometric, approach. In Ref. |44|, Coniglio showed that the p-derivative of Pp,G[2; y] for two 
vertices x,y G V (not necessarily neighbors) satisfies the following differential equation: 

Epp G[a: GGy] = ^Ep^G[K,y], (4.1) 

where Xx^y{A) is the number of bridges on anjQ self-avoiding path in A connecting x and y. We emphasize 
that for a given edge (x, y) G E the quantity \x,y can be non-zero, even if (x, y) ^ B{A). This is because 


^It can be verified that the number of bridges on a self-avoiding path between to vertices is an invariant, hence it does 
not matter which self-avoiding path one chooses (in case more than one exists). 
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Xx^y counts the number of bridges in the spanning subgraph (V, A) that would disconnect x from y upon 
removal. As an example, consider vertices b and e in the left panel in Figure Even with (6, e) not 
being part of the considered spanning subgraph, we have Xb^e = 2, because edges (in fact bridges) (6, c) 
and (c, /) are essential for the connectivity of b and e. Note that bridge (a, b) is not essential for the 
connectivity of b and e as its removal does not disconnect b and e. Note that relationship (4.11 also 


follows straightforwardly from the Russo-Margulis formalism applied to the indicator function of the 


event x y. As a result of the bijection relation (3.11 we can transform (4.1) into a differential equation 
for the bridge density: 






(x,y)eE 


1 

p{l-p) 


^p,g[B]. 


(4.2) 


We emphasize that now Xx^y 
{x,y) G E. Observe 


is only evaluated for nearest-neighbor pairs, that is for x,y € V with 




{x,y)^E e^E 


A 




= E 

eeE 


E *{ 

x,y)eE '■ 


A 

x-^y 




(4.3) 


Now, for a given edge e G E define PeiA) as the sum in the square brackets in (4.3) or, equivalently, to be 
0 when e is not a bridge in A and otherwise the number of nearest-neighbor pairs for which the bridge e 
lies on any (all) self-avoiding path(s) in {V,A) between them. In other words PeiA) counts for how many 
nearest-neighbor pairs edge e contributes to Xx,y We call pe the bridge load of e. Note that we have 
Pe(A) > l{eGB(A)} and hence Ep_q[yOe] > f'p,q[e G B], For an example consider the left panel in Figure]^ 
here edges (a, 6), (6, c), (c,/) have bridge loads 2, 3, 3, respectively. It is evident that in a certain sense 
Pe generalizes the bridge density and encodes how relevant a given bridge is for the connectivity at large. 
We emphasize that one only has to count connected neighbor pairs. We finally obtain: 




1 


p{l-p) 


]Ep,G [B], 


(4.4) 


where we defined p{A) = EeeE /Oe(A). Therefore we find that the l.h.s. in (4.4) vanishes for p = pf, 
where pf is the solution of the following equation 


^P{,g[B] = PfE. 


'Pi,G 


[P\- 


(4.5) 


We remark that pi is a graph dependent quantity and we suppress the explicit dependence. In Fig. 
we illustrate this condition for a small system size and the case G = Moreover, under translational 
invariance (which is given, e.g., for Z|^), (4.5) is equivalent to Ppj._G[e G B] = pfEp^^Q[pe] for arbitrary 
e G E. This has a nice intuitive interpretation: Ep^ciPe] counts the typical number of nearest-neighbor 
pairs at the interface between two clusters that are glued together only by the presence of bridge e. There 
are two effects at work here: increasing p typically increases the size of the clusters and therefore possibly 
also the number of nearest neighbors on the interface. Additionally, however, increasing p will eventually 
create additional links between the two clusters, in which case e ceases to be a bridge, such that this 
effect decreases Pp,G[e G B]. Now, the maximal fragility (maximum of Pp,G[e G B]) is attained for the 
bond density pf, where the expected number of closed edges between two clusters glued together by e 
equals precisely the probability that e is a bridge. Increasing p further strengthens the cohesion between 
the clusters, and hence the probability of e being a bridge starts to decrease afterwards. 

In what follows we restrict the analysis to the case G = in order to extract the asymptotic value of 
Pf for Z^. Consider Figure]^ where we numerically confirm condition (4.5) and estimate pf = 0.4056(5) 
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Figure 6: Numerically estimated density of bridges and bridge load for percolation on for L = 16 and L = 128, 
corresponding to the smaller and larger markers, respectively. Additionally, the inset shows our estimates for pf for system 
sizes from L = 8 up to L = 128. The estimates for the location of the crossing “saturate” to a value 0.4056(5), and therefore 
we conclude that pf = 0.4056(5) for Z^, according to ( |4.5| |. We point out that we used a naive depth-first search traversal 
based algorithm for determining the bridge load, so we could not study systems with L > 128 in practise. The lines in the 
figure are obtained by interpolating the data points using splines. 


for We note that, albeit the maximum of Pp,z 2 [e £ B] is attained for p = pf, the bridge load [pe] 
increases further and is apparently maximized for p = Psd(l) = Pc — 1/2- This nicely reflects the intricate 
cluster structure at criticality. One might wonder how the bridge density can decrease but the overlap 
continues to increase. This can happen when in most (probable) configurations e is not a bridge, but 
for configurations where e is a bridge one has typically a very large self-entanglement of the clusters 
attached to the two ends of e, outweighing the effective decrease of Pp z 2 [e G B], This effect is most 
drastic for p = p^. We note that Coulomb gas arguments predict that the l.h.s. of (4.4) remains finite 
at p = Psd(l) = 1/2, which clearly implies that Ep ^^iPe] remains finite at p = 1/2. More precisely we 
can show that the peak of E^ 2 |[Pe] scales as ^ 2 - 22:2 £pj. q g [0,4], where in particular X 2 = 5/4 for 
percolation [see the discussion in the next section and in particular (5.3)]. Unfortunately, so far we have 
not been able to establish an exact value for pf(l) let alone in the correlated setting 5 ^ 1 for Z^. This 
is a challenging problem we wish to address in a forthcoming study. 


5. Bridge fluctuations 


As we have seen, the bridge-edge formula ( 2 . 12 ) together with the theory of finite-size scaling, provides 


a rather complete understanding of the critical bridge density. Higher moments of the bridge distribution 
can be discussed with similar techniques as we will show now for the example of the variance. Naively, 
one might expect that this variance Varp^q^c]]!?]] only depends on the fluctuations of N. In this case the 
critical variance would follow Var^ ^[jHj]/L‘^ « which, in two dimensions, in turn would imply a 

divergence with L for g > 2 and a “saturation” for g < 2 [45]. As we will show however, the story is not 
quite as simple. 


To work out the fluctuations of the bridge density, we apply the Russo-Margulis formula to the second 
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derivative of the partition function Zp g Q = q, G), which we then in turn equate to the expression 

one obtains by explicit differentiation. To start with, we have, using (2.11), 


dlZp^g^G = (1 - q)dp 


= -Zp,g,G^^p,,,G [|S|] + ^dpEp^G [q^\B\] . 


Let us now focus on the second term: 

^dpEp^G [q^\B\] = — E [^e {q^MfGB})] 


eJ^E 


1-9 




(5.1) 


eJ&E ace: 

eGA 


We can split the inner sum into two sums corresponding to e G B{A) and e € G(A), of which we first 
consider the former 


1-9 


X! '^P,g[M l{e=/}9^^'^^l{/GB(A)} + l{e#/}9^^^Hl - 9)1{/GB(A)} 


ejGB 

c<=B(,A) 


A few comments are in order: If e = / we clearly have that e = / ^ B{Ae) hence we recover the 
expected bridge density. For e ^ / it is important to observe that removing a bridge e cannot influence 
the pivotality of an occupied edge /, in particular when / G B(A) then also / S B(Ag), whenever e is a 
bridge in A. The above can therefore be re-written as 


1-9 


Zp,g,G^p,g,G [|1?|] + - -^ Pp. 9 ,G[e G B,f G B], 


e^fdE 


For the non-bridge sum of (5.1) note that only summands with e ^ / contribute, because otherwise if 
e G C(A) then trivially / = e ^ B{A), B{Ae). Furthermore, if e ^ / such that e G C(A) and / G B(A) 
then we have also a vanishing contribution because deleting a non-bridge cannot change the fact that / is 
a bridge. Thus, there can only be a contribution for a configuration A such that e G C(A) and / G C(A) 
as well as both edges are in one cycle and deleting e will destroy the cycle and hence cast / into a bridge 
in Af,. Moreover, both edges e and / must only be in one cycle that has no edge overlap with other cycles 
(imagine two clusters glued together in parallel by the edges e and /, hence in particular there are no 
additional links between the two clusters). Write e <4> / for the above event involving edge e and /. We 
obtain for the second term: 

~Zp^g^G y ' IPp.ij.G [g /] . 




e^/GE 


We therefore obtain eventually for d^Zp^g-. 


dpZp^q^G = Zp^g^G —^ ( (1 - 9 ) X! '^P.9.G[e & B,f G B]- ^ lPp.<j.G[e /] 
P ' e^feE e^feE 


= z 


1-9 


P:<1:G 2 


(1 - q) Varp,,,G[|-B|] -b (1 - q)Ep^g^G[\B\]‘^ - (1 - 9 )Ep, 9 ,G[|^|] - ^P.i.cie f] 

e^fGE 
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On the other hand one can show by explicit differentiation and after some straightforward but tedious 
algebra that one has: 


1 Ep,g,G[N]{2p-l-2p\E\)-\E\p'^ + \E\^p'^+Vaip,g,G[N]+Ep,g,G[N]^ 


-‘p,q,G 


-dlZ, 


p,q,G 


p 2 (l — pY 


Equating both expressions for diZp^q^G with subsequent rearranging yields 


Yaip._q^G[\B\] = 


Ep,,.G[iV] (2p - 1 - 2p|£;|) - \E\p'^ + Varp,,,G[iV] + (Ep,,.G[Af] - \E\pf + 2|£;|pEp,,,G[iV] 


{l-qYd-pf 


1 


]Ep,5,g[|^|] + ]Ep,(J,g[|^|] + , _ ^p,q,GY ^ f]- 

® e^feE 


Using the bridge edge identity (2.12) we can simplify the above to 

Varp_q_G[|-S|] _ Ep,,,g[A/ 1 (2p - 1) - + Varp_q,G[^]/l-E| 


\E\ 


(1 - g)2(l -p)2 


+ Ep,g,G[^] + _ 1^1 IPp.g.G[e ^ /]• 


e^fGE 


(5.2) 

We emphasize that (5.2) is an exact result valid for any p and q as well as any graph G, thus it has the 
same status as the bridge-edge identity (2.12). Furthermore, we remark that for e = (it, v) and / = (x, y) 
we have by a bijection argument, similar to the one used in Sec. [^to derive the relationship between the 
bridge density and nearest-neighbor connectivity: 


d,9,g[6 ^ /] — 


p^ 1 


(1 -p )2 q 


2 „ 


,g[(m,'c) ^ {x,y)]. 


Here Pp_g^G[(w,n) M (x,y)] for (it, u), (cc, p) S if is the probability that the two nearest-neighbor pairs 
belong to two different clusters, such that the two distinct clusters each contain one vertex from {it, i;} 
and one vertex from {x,y} (see Fig. 1(c) in Ref. [48]). 

We pause for a moment, and remark that completely analogous arguments can be used to obtain an 
expression for the p-derivative of Ep_q_G[^]) studied in Section]^ in terms of the probabilities of events 
e /. A consequence of this is that we can obtain an explicit expression for Ep_q_G[p]: 

Ep,(3,g[p] = Ep,q,G[^] + ~ IPp,q,G[e /]. (5-3) 

' ' ^ e^f&E 


In particular, this shows explicitly that the bridge load can increase further beyond the point pi{G), for 
which Ep^g[B\ is maximal. Moreover, as we show below, we find that at criticality, in two dimensions, 
[p] ~ which together with X 2 = 5/4 for the case <7 = 1 in two dimensions, shows that the 

peak at p = 1/2 in Figure]^ stays bounded as L —>■ 00 . 

Coming back to the study of fluctuations of \B\, we focus in what follows on the self-dual line for the 
RC model on and study the continuum limit. As shown by Vasseur et al. |48j one has for two pairs 
of neighboring vertices (it, v) and (x, y) at distance r the following asymptotic for large r 

^q,i?[{u,v) y {x,y)\ ~ A( 9 )r"^“="(«). 


This follows from the construction of a four-leg watermelon event, due to the four hulls propagating 
from the neighborhood of {u,v) to the neighborhood of {x,y), which are associated to the two clusters 
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involved. We hence expect the following asymptotic behavior 


H ^ H [(w, v) y (x, y)] 

' ' ey/6B I I iu,v)^{x,y)eE 

L 

« A(q) y ' dr27rrr-2^=^(«) 

= a(g)L 2 - 2 x.(,)+^(^)^ (5.4) 


where a{q) and /3(g) are two q dependent constants. 


Inspecting the form (5.21, we hence see that the fluctuation of the bridge density has two contributions, 


one related to the variance of the edge density that scales as at criticality, and another one related 
to the above mentioned watermelon event with scaling proportional to ^ 2 - 22:2 order to decide which 
contribution constitutes the dominant leading scaling, it is useful to recall some relevant exact results from 
the Coulomb gas approach to two dimensional critical phenomena [1]. One finds that ajv = 4 — 12/g and 
a ;2 = 2 — (4 —g)(4 + 3g)/8g, where for critical 2D RC models one has to relate the Coulomb gas coupling g 
to the cluster weight q via g = 2 + 2 cos (g 7 r/ 2 ) [3^. A small calculation shows then that constitutes 
the leading contribution to the system size scaling of the variance whenever g > 1. In particular, for 
g > 2, where ajv > 0, this translates into a divergence of Var^ z 2 [|i3|]/L^ with L, resembling the specific- 
heat singularity in the same cluster weight regime. On the other hand, for g < 1 the leading term 
is proportional to Interestingly (and somewhat unexpected), one can explicitly verify that for 

q < q = 4cos^ ( 7 r/-\/ 3 ) = 0.2315891 •• • the variance of |i3| diverges with the exponent 2 — 2 x 2 > 0. 
This particular value g follows from the condition 2 — 2 x 2 (g) = 0 X 2 (g) = 1, which, using the above 

expressions, yields a quadratic equation for g = g{q) (or g): 8 g = (4 — g)(4 -|- 3g), with the two solutions 
g = ±4/-\/3. Taking only the positive solution we find that g = 2 -|- 2 cos {2^T/^/i) = 4cos (tt/i/S) . The 
situation is summarized in the dotted and solid lines of Fig. showing the Coulomb gas values of ajv 
and 2 — 2 x 2 , respectively. We also analyzed the variance of the bridge density numerically. The fitting 
functions and resulting parameters are summarized in Table[^ and the corresponding parameter estimates 
are indicated by the symbols in Fig. Clearly, we find excellent agreement with the expectations from 
Eq. (5.2) discussed above. For the marginal value g = 0.2315891 • • • we expect a logarithmic divergence. 


and we indeed find the corresponding form to yield the best fit to our simulation data. 


It remains to discuss the percolation case g = 1 where the bridge-edge identity (2.121 becomes singular 


and the above derivation needs to be revisited. Here, we derive the singular behavior based on another 
bijection argument that allows us to harness recent results on logarithmic observables emerging from a 
careful analysis of the appropriate logarithmic conformal field theory description of critical percolation 
[23148]. Before we do so, note that in order to extract the asymptotic scaling of the variance it suffices 
to study the (co-)variance = IPp.g.G[e S H,/ G H] — Pp,q,G[e G B]Fp^q^c[f G B], which relates to 
Varp_q_G[|i3|] via the well known identity 


Varp,,,G[|B|] = qej- 

eJeE 

To start with, recall that we consider critical bond percolation on Z|, i.e. g = 1 andp =1/2 and write P[-] 
for Pi/ 2 ,i,z| [•]• Furthermore, note that is a transitive graph, and hence none of the following events 
depend on the explicit edge or vertex, used in the arguments. Now, fix two edges e = (xi, j/i), / = (x 2 , 2 / 2 ) 
that are distance r <C T apart. Consider the event {e G H A / G B}. All configurations contributing 
to this event can be further sub-divided into two events, depending on whether e and / belong to the 
same connected component in {V,A) or not. Denote the two events by Hi and D 2 , respectively. Choose 
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Figure 7: The exponents a/u and 2 — 2x2 appearing in the system size scaling of Varg[|S|]/m at criticality. The solid and 
dotted lines show the exact value of 2 — 2x2 and a/u, respectively, following from the Coulomb gas mapping. The symbols 
denote our numerical estimates from fitting to the variance of the bridge density, cf. the data collected in Table 


a configuration A that belongs to ^ 2 , i-e. the two edges e and / are bridges in (V, A) and belong to two 
different connected components. The crucial point is that we can relate A one-to-one to A — {e,/}, a 
configuration where Xi, yi, X 2 , 2/2 belong to four different clusters. Denote all configurations in which the 
four vertices belong to four different components by CI 2 (this is a event). Due to the choice of g = 1 and 


p = 1/2 we have that P[n 2 ] = ■ 


D2 


Let us now consider the event fli, i.e. the set of configurations for 


which e and / belong to the same component and both e and / are bridges. Now, because both edges 
are pivotal we can relate any such configuration A £ VLi one-to-one to a configuration where Xi and yi as 
well as X 2 and 1/2 are disconnected, and for which xi,yi,X 2 ,y 2 belong to three different clusters, of which 
one cluster contains one vertex of {xi, yi\ and one of {x 2 ,y 2 \- Denote the corresponding event fli. Note 
that any configuration A G where e and / are bridges belonging to the same cluster, must yield 3 
disconnected clusters in A' = A — {e, /}. This is because the alternative case of 2 disconnected clusters 


P[f2i] = 


111 


Finally, the probabilities 


3 in 

A, which is 

1 

and P 

O 2 


T 

lun 


were studiec® in 05] in the framework of a 


corresponding logarithmic conformal field theory. We note that in the continuum limit both probabilities 
P[ni] and P[n 2 ] only depend, to leading order, on r. Since we have P[e G B,f G B] = P[f2i] -|- P[n 2 ] we 
obtain by falling back to [48] : 


r]e,f ~ [o-I-6 log (r)]r 


— 2X2 


Vari/2,i,z| [l-S|] 

Z2 


+ [ 6 ' + c'log (L)] 


(5.5) 


where X 2 = 5/4 is the two-arm exponent for critical percolation and a, 5, a', b' and c' are constants. 
Our numerical analysis confirmed this scaling, as shown in Table and Figure Thus considering the 
variance of bridges for critical percolation yields yet another manifestation of the underlying logarithmic 
conformal field theory [48] . We remark that the same logarithmic multiplicative corrections are expected 


®Vasseur et al. in [48] denote P[ni] by Pi(r) and P[f22] by Po(r). 
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q 

-^min 

xVN^.o.i. 

Model 


afv 

2 — 2x2 

0.0005 

12 

1.01 

a + hL'^ 

0.462(16) 

-1.9576 

0.4752 

0.005 

4 

1.15 

a + hL^ 

0.422(3) 

-1.8679 

0.4222 

0.05 

16 

0.58 

a + hL'^ 

0.258(5) 

-1.6005 

0.2599 

0.1 

8 

0.76 

a + hL^ 

0.164(4) 

-1.4492 

0.1648 

0.231589 

8 

1.01 

a -1- 6 log (L) 

- 

-1.1962 

0 

0.5 

12 

1.4 

a + bL^ + 

-0.217(51) 

-0.8778 

-0.2191 

1 

6 

0.7 

a + bL'^ [1 -1- clog (L)] 

-0.496(21) 

-0.5 

-0.5 

1.5 

4 

1.05 

a + -f cL 2 - 2^2 

-0.218(15) 

-0.2266 

-0.7205 

2.5 

48 

0.94 

a-f 

0.202(14) 

0.2036 

-1.1052 

3 

6 

1.14 

a-f 

0.398(3) 

0.4 

-1.3 


Table 2: Numerical results for the leading exponent k, in the finite-size scaling of Var^ ^2 [|-B|]. The two rightmost columns 
show the exact values obtained from the Coulomb gas mapping. The exponent values in bold face indicate the asymptotically 
dominant behavior. For the cluster weights q = 1.5 and q = 0.5 we also performed a fit to the form a + bL'^^ which yielded 
slightly worse results, due to the proximity of the two (negative) exponents. In all cases the quality-of-fit Q was at least 
5%. 


in higher dimensions [49) . More precisely, one expects even the same scaling form, that is 

Varp^^i y[|g|] ^ ^ ^ 


where we used the well known fact that X 2 = d — 1/v holds for critical percolation [44]. It would be 
interesting to to verify this theoretical prediction numerically. 

Before we proceed to the next section we briefly discuss the finite-size scaling of Var^^^i [|C|]/|i?|. This 
can be worked out by observing that \C\ = |^| — \B\ which in turn implies the general identity 

^&'^p,q,G[\C\] = Varp,,,G[|^|] + Varp,g,G[|-S|] - 2CoVp,,,G[|^|, |^|]- 

Additionally, one can verify that (see m or by explicit differentiation of Zp^q) 

Covp,,,G[|^|, 1^1] = pil - p)dpEp^q_G[\B\]. 

An immediate consequence of these observations is that Var^ ^hl^H governed by the same finite-size 
scaling as Va.iq j^d[\B\] and we therefore omit further details. 


6. Bridge-free clusters 

In concluding the present study, we return to the discussion of the structure of the incipient percolating 
cluster outlined above in the introduction. In the ‘links-nodes-blobs’ picture of Ref. [T5| the attention 
is focused on the backbone of the cluster as the skeleton ensuring long-range connectivity. Bridges on 
the backbone are singled out as the bonds carrying full current for the case of an external potential 
difference being applied, hence their name ’red bonds’. Due to the importance of these structures for 
a number of applications of percolation theory, significant effort has been invested in the determination 
of the associated fractal dimensions deB and (Irb of the backbone and the red bonds, respectively [50]. 
Even in two dimensions, where most critical exponents and fractal dimensions, including dp and d^Bj are 
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Figure 8: Effective backbone scaling exponents according to Eq. ( |6.1[ | for the q = 1 RC model extracted from the scaling 
of the backbone sites, the backbone bonds (data from Ref. [50]), the junction-free clusters and the bridge-free clusters (our 
data), respectively. 


known exactly from Coulomb gas and further arguments [3], there is no exact expression for the fractal 
dimension of the backbone. 


Clearly, the notions of backbone and red bonds appear to depend on the concept of long-range 
conductivity, and it is not obvious how they could be defined in terms of the more local connectivity 
structure of bridges and non-bridges. The concepts of junctions and bridges introduced in Ref. [TH] are 
related, but not commensurate to the idea of separating bridges in the backbone (red bonds) and bridges 
in the blobs. We expect that, away from g —>■ 0, most bridges in the backbone are junctions as the ‘links- 
nodes-blobs’ picture tells us that the backbone has cycles. Still, there are junctions in the dangling 
ends and there are branches in the backbone. Also, as Xu et al. established for percolation in Ref. [T5] . 
branches and junctions are asymptotically finite fractions of the bridge set and hence of the edge set. It 
is therefore interesting to study the fractal dimensions of the related edge classes. Removal of branches 
from critical configurations will typically only shave off the last part of dangling ends, and we do not 
expect this to alter the fractal dimension of the incipient percolating cluster. This is indeed what Xu et 
al. report [18] and we find the same to be true also for the RC model. On the other hand, removing 
the junctions means deleting most of the red bonds and hence breaks down the percolating cluster into 
individual blobs. Their size scales with the backbone fractal dimension as oc . The same holds true 
for the case of bridge-free clusters where both, branches and junctions, have been removed. Studying 
these sets hence provides an alternative route to the determination of cIbb- 

To investigate the utility of this approach for the determination of cZbb, we determined the scaling 
of the junction-free and bridge-free clusters for the percolation case q = 1. If we consider the effective, 
system-size dependent fractal dimensions m 


dsB.efi 


\og[N{2L)/N{L/2)] 
log 4 


( 6 . 1 ) 


where N{L) denotes the number of sites in the junction-free or bridge-free clusters at size L, respectively, 
we can easily compare corrections to scaling for the different approaches. This is illustrated in Fig. 
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Table 3: Fit results for backbone fractal dimension (Iqb. 


q 


-^min 

A^d.o.f. 

X^/^d.o.f. 

Q 

0.0005 

1.264(2) 

32 

8 

1.0631 

0.3858 

0.0050 

1.301(1) 

16 

10 

0.2934 

0.9829 

0.0500 

1.383(1) 

16 

10 

0.4661 

0.9126 

0.1000 

1.425(1) 

12 

11 

1.1773 

0.2966 

0.2000 

1.480(1) 

12 

11 

1.0047 

0.4390 

0.3000 

1.516(1) 

32 

8 

0.7543 

0.6434 

0.5000 

1.567(1) 

24 

9 

0.9153 

0.5104 

0.7000 

1.603(1) 

32 

8 

0.6676 

0.7206 

0.9000 

1.630(1) 

24 

9 

1.0700 

0.3813 

1.0000 

1.643(1) 

24 

9 

1.3227 

0.2187 

1.2500 

1.670(1) 

24 

9 

1.0177 

0.4227 

1.5000 

1.692(1) 

24 

9 

1.1271 

0.3389 

1.7500 

1.716(2) 

128 

4 

0.7496 

0.5581 

2.0000 

1.732(1) 

96 

5 

1.2095 

0.3016 

2.2500 

1.744(1) 

48 

7 

1.0578 

0.3880 

2.5000 

1.756(1) 

96 

5 

0.5469 

0.7408 

2.7500 

1.774(1) 

96 

5 

1.1302 

0.3417 

3.0000 

1.785(1) 

64 

6 

0.6563 

0.6851 

3.2500 

1.797(1) 

96 

5 

0.9530 

0.4452 

3.5000 

1.806(1) 

96 

5 

0.8288 

0.5289 

3.7500 

1.828(4) 

400 

2 

1.0609 

0.3462 

4.0000 

1.83(1) 

32 

7 

1.3921 

0.2035 


comparing the exponents extracted from the scaling of the junction-free and bridge-free clusters to the 
effective exponents found from the scaling of the backbone itself [50]. Clearly, the scaling of bridge- 
free clusters is significantly less affected by scaling corrections than the scaling of the backbone itself. 
Compared to this, junction-free clusters show slightly increased corrections, however interestingly with a 
correction amplitude of the opposite sign. We remark that our estimate for percolation cIbb = 1.643(1) is 
consistent with existing literature values, such as dee = 1.64336(1), 1.6434(2), 1.6432(8) in [TTlIT^ 1^ . 
respectively. As this is a small-scale study not using the particularly efficient algorithms available for the 
uncorrelated percolation case, the reduced scaling corrections do not immediately lead to an improved 
estimate. Instead, we focus on the estimation of dsB for a wide range of q values for which no estimates 
have been previously reported m- Whereas for the determination of des directly from the backbone 
corrections to scaling need to be carefully taken into account, it turns out that the bridge-free clusters 
provide a scaling route less encumbered with such problems. As is evident from the results collected in 
Tablej^ over the full range of q values we achieve excellent fits without the inclusion of scaling corrections. 
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Figure 9: Estimated backbone fractal dimension using the method described in m- For comparison also values from na 
are shown. The disagreement for g = 4 is most likely due to known logarithmic corrections. 


7. Discussion and outlook 

As we have seen, the natural classification of edges into bridges and non-bridges provides an under¬ 
standing of how the cluster weight influences the connectivity in the random-cluster model. For any 
g ^ 1 we established that both bridges and non-bridges have expected densities that are linearly related 
to the overall density of open edges, such that there is always a non-zero fraction of both edge types. 
The derivation is based on an application of the Russo-Margulis formula to the partition function of the 
model, a connection that allows to express analytical derivatives in terms of combinatorial and geomet¬ 
ric quantities. For the percolation limit g —> 1 the linear relation does not hold and rather transforms 
into an intuitive identity relating the bridge density to the covariance between the number of connected 
components and open edges. The latter identity explains previously observed finite-size corrections for 
this covariance in the percolation model |35]. Thus, the established bridge-edge formula allows us to 
connect the finite-size corrections of the densities of bridges and non-bridges to corrections in the energy 
density governed by the exponent d — Ijv. The lack of finite-size effects in the density of the so-called 
type-1 edges for the special case of a 2D lattice with periodic boundary conditions, corresponding to a 
torus graph, can be understood from a cancellation of correction terms, and we clarified the character 
and origin of correction terms with exponents d — Ijv and —X 2 ^ respectively [18) . 

By studying the dependence of the bridge density on the bond occupation p we observe that its 
maximum is attained for pi{q) strictly below the corresponding critical point Psd(g)- We explained this 
observation for the percolation case by studying the newly introduced bridge load of an edge. For 
percolation on we estimate pi = 0.4056(5). Our numerical study of this quantity was limited by the 
computational cost of determining the bridge load for each edge, which we were not able to calculate 
with linear (worst-case) computational effort for all edges. However, we remark that it is possible to 
express a condition for pf in terms of easily computable covariances and moments (in linear time), which 
however requires to estimate the difference of quantities of order 0{L‘^) to a precision of order 0(1), which 
implies a statistical inefficiency. The construction of a more efficient algorithm for the determination of 
the bridge load constitutes an interesting open problem. 
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Having established the expectation values of the bridge and non-bridge densities, we turned to a 
study of the fluctuations. Again using the Russo-Margulis formalism we established a second bridge-edge 
identity relating the variances Varp_q_G[|i?|] and Varp_g_( 3 [|C'|] to the variance of the bond density and an 
extra term relating to a specific four-leg watermelon event [48] . For the special case of the two-dimensional 
model, we use results on the scaling of this event derived from conformal field theory [48j to predict a 
singularity in both variances for q < q = 4cos^ (7r/-\/3) = 0.2315891 • • •, an effect completely absent in 
the fluctuations of the overall number of open edges. Numerical simulations confirm these predictions 
to a high precision. The limiting percolation case required a separate analysis, and here we are able to 
show that the finite-size corrections to the density of bridges for critical percolation have logarithmic 
multiplicative corrections, a direct consequence of the underlying logarithmic conformal field theory. It 
would be interesting to study these fluctuations for higher dimensional systems to see whether such a 
singularity appears generically for some q < 1, and how this threshold value depends on dimensionality. 
It follows from hyperscaling, vd = 2 — a^ and the relation l/u = dn = d—X 2 |H)] valid only for percolation 
(g = 1) that in this case ajv = d — 2x2, irrespective of dimension d. Hence, g = 1 is a marginal case 
in all dimensions below the upper critical one, yielding a multiplicative logarithm in the variance of the 
number of bridges. It would be interesting to see whether d — 2x2 > otjv for g < 1 in 3D also (see 
[45] and [45] for two definitions of X 2 that can be generalized to d > 3), thus hinting at an interesting 
physical/geometrical difference between the (anti-monotone) g < 1 and (monotone) g > 1 regimes. 

Removing bridges leads to fragmentation of clusters as discussed recently in Ref. [20] • Removing all 
bridges decomposes the percolating cluster into blobs, and it turns out that this approach allows for a 
very precise determination of the backbone fractal dimension, much less affected by finite-size corrections 
than the more traditional approach of studying the backbone directly [50] . We provide estimates of this 
important dimension, which is one of the few exponents of the random-cluster model in two dimensions 
for which no exact expression is known, for a wide range of values in 0 < g < 4, many of which have been 
studied here for the first time. 

An even more detailed understanding of the role of bridges could be expected from a generalization 
of the random-cluster model giving different weights to bridge- and non-bridge bonds. This problem 
of a generalized bridge percolation [TH] is the subject of a forthcoming study. The 0{n) loop model is 
another general class of frequently studied systems m- As these are graph polynomials as well (namely 
over a restricted set of spanning subgraphs, i.e., all loop configurations or Eulerian subgraphs [52]), it is 
therefore tempting to investigate in how far a similar Russo-Margulis approach can be applied there or, 
more generally, to Eulerian subgraph models [52] . possibly yielding novel identities. 
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Appendices 

A. Influence of an edge 


In order to prove (2.8), we first decompose the sum over subgraphs AGE into two contributions. 


corresponding to configurations for which e € A and the complementary set of configurations with e ^ A. 
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Then observe that for A such that e € A we clearly have A^^ = A and similarly for A with e ^ Awe have 
A^ = A. Furthermore we can also relate the probabilities of A and the modified set Ag or A^ once we 
know whether e G A. Using these rules, we arrive at 


]Ep,G[<5e-^] 


^ Pp.g[A]{X(A^)-X(A,)} 

ACE 

Vj,,g[A]{X{A)-X{A,)}+ Y, Vj,,g[A]{X{A^) - X{A,)} 

ACE: ACE: 

eGA e^A 

Y Pp.g[A]{X(A)-X(A,)} + 1^ ^ F,^a[A]{X{A)-X{A,)}) 

ACE: ^ ACE: 

eGA eGA 

^ ^ F,^g[A]{X{A)-X{A,)}. 

P ACE: 
eCA 


We remark that in the last sum one can safely discard the condition e S A, because for A C if with e ^ A 
one trivially has Ag = A and hence X{A) — X{Ae) = 0. 


B. Free energy derivative 


In order to evaluate Ep_G[9^], note that K only changes on removing a bond e if it is a bridge, i.e., 
K{A) — K{Ae) = —l{eGB(A)}- Hence, we have 

qK{A) _ qK{A^) ^ qK{A)f^^ _ 


We can now apply (2.8 1 to obtain: 

F^,G[5eq^\ = - E ^P.g[A] 


Hence we obtain 


qK(A) _ ^K(A^) 


ACE: 

eGA 


l-q 


Y ^P,G[A]q^^^"''^{eCB{A)} 


ACE: 

eGA 


l-q 

P 

l-q 

P 


Zrc{P, q, G')Ep_g_G[l{eGB}] 
Ep.G[<?^]Pp.9.G[e e B]. 


_5 

dp 


— \ogZnc{p,q,G) = — logEp,Gb ] 


_5 

dp 


1 d 


^P,Gk^]dp ^ 


Ep,G[g 


K-i 


1 


1-9 


E ®P:G[<5e 


Wl 


P 


E/ ^P>9,G[e G B] 


eCE 


I^^I^Ep.,,g[H]. 


28 













References 


[1] D. Stauffer, A. Aharony, Introduction to Percolation Theory, 2nd Edition, Taylor &: Francis, London, 1994. 

[2] G. Grimmett, Percolation, 2nd Edition, Springer, Berlin, 1999. 

[3] B. Bollobas, O. Riordan, Percolation, Cambridge University Press, Cambridge, 2006. 

[4] B. Nienhuis, Coulomb gas formulation of two-dimensional phase transitions, in: C. Domb, J. L. Lebowitz (Eds.), Phase 
Transitions and Critical Phenomena, Vol. 11, Academic Press, London, 1987, p. 1. 

[5] J. L. Cardy, Conformal invariance, in: C. Domb, J. L. Lebowitz (Eds.), Phase Transitions and Critical Phenomena, 
Vol. 11, Academic Press, London, 1987, p. 55. 

[6] O. Schramm, Scaling limits of loop-erased random walks and uniform spanning trees, Israel J. Math. 118 (2000) 221. 

[7] S. Smirnov, W. Werner, Critical exponents for two-dimensional percolation. Math. Res. Lett. 8 (6) (2001) 729—744. 
doi:10.4310/MRL.2001.v8.n6.a4 

[8] S. Smirnov, Critical percolation in the plane: Conformal invariance, Cardy’s formula, scaling limits, C. R. Acad. Sci. 
Paris Ser. I Math. 333 (3) (2001) 239-244. 

[9] H. Duminil-Copin, S. Smirnov, Conformal invariance of lattice models, in: D. Ellwood, C. Newman, V. Sidoravicius, 
W. Werner (Eds.), Probability and Statistical Physics in Two and More Dimensions, Vol. 15 of Clay Mathematics 
Proceedings, American Mathematical Society, Providence, RI, 2012, pp. 213—276. 

[10] C. M. Fortuin, P. W. Kasteleyn, On the random-cluster model i. introduction and relation to other models, Physica 
57 (4) (1972) 536-564. doi:10.1016/0031-8914(72)90045-6 

[11] G. Grimmett, The random-cluster model. Springer, Berlin, 2006. 

[12] F. Y. Wu, The Potts model, Rev. Mod. Phys. 54 (1982) 235-268. doi: 10.1103/RevModPhys. 54.235 

[13] V. Beffara, H. Duminil-Copin, The self-dual point of the two-dimensional random-cluster model is critical for g > 1, 
Probab. Theory Rel. 153 (3-4) (2012) 511-542. doi: 10.1007/s00440-011-0353-8 

[14] R. J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press, London, 1982. 

[15] H. E. Stanley, Cluster shapes at the percolation threshold: and effective cluster dimensionality and its connection with 
critical-point exponents, J. Phys. A 10 (11) (1977) L211. doi: 10.1088/0305-4470/10/11/008 

[16] A. Coniglio, Fractal structure of ising and potts clusters - exact results, Phys. Rev. Lett. 62 (1989) 3054. doi: 
10.1103/PhysRevLett.62.3054 

[17] Y. J. Deng, H. W. J. Blote, B. Nienhuis, Backbone exponents of the two-dimensional g-state Potts model: A Monte 
Carlo investigation, Phys. Rev. E 69 (2004) 026114. 

[18] X. Xu, J. Wang, Z. Zhou, T. M. Garoni, Y. Deng, Geometric structure of percolation clusters, Phys. Rev. E 89 (1) 
(2014) 012120. 

[19] N. A. M. Araujo, K. J. Schrenk, J. S. Andrade, H. J. Herrmann, Bridge percolation, Preprint arXiv:1103.3256 arXiv: 
1103.3256 

[20] E. M. Elgi, M. Weigel, N. G. Fytas, Fragmentation of fractal random structures, Phys. Rev. Lett. 114 (11) (2015) 
115701. 

[21] R. H. Swendsen, J. S. Wang, Nonuniversal critical dynamics in Monte Carlo simulations, Phys. Rev. Lett. 58 (1987) 
86-88. doi:10.1103/PhysRevLett.58.86 

[22] L. Chayes, J. Machta, Graphical representations and cluster algorithms II, Physica A 254 (1998) 477. 

[23] E. M. Elgi, M. Weigel, Efficient simulation of the random-cluster model, Phys. Rev. E 88 (2013) 033303. doi: 
10.1103/PhysRevE.88.033303 

[24] M. Sweeny, Monte Carlo study of weigted percolation clusters relevant to the potts model, Phys. Rev. B 27 (1983) 
4445. 

[25] R. G. Edwards, A. D. Sokal, Generalization of the Fortuin-Kasteleyn-Swendsen-Wang representation and monte carlo 
algorithm, Phys. Rev. D 38 (1988) 2009. doi: 10.1103/PhysRevD. 38.2009 

[26] J. Salas, A. D. Sokal, Dynamic critical behavior of the Swendsen-Wang algorithm: The two-dimensional three-state 
Potts model revisited, J. Stat. Phys. 87 (1-2) (1997) 1-36. doi: 10.1007/BF02181478 

[27] H. Hu, H. W. J. Blote, R. M. Ziff, Y. Deng, Short-range correlations in percolation at criticality, Phys. Rev. E 90 (4) 
(2014) 042106. 

[28] G. Grimmett, Probability on Graphs - Stochastic Processes on Graphs and Lattices, Cambridge University Press, 
Cambridge, 2010. 

[29] J. E. Steif, A mini course on percolation theory, Jyvaskyla Lectures in Mathematics 3 (2011) 1. 

[30] J. W. Essam, Connectedness and connectivity in percolation theory, in: A. Barlotti, M. Biliotti, A. Cossu, G. Korch- 
maros, G. Tallini (Eds.), Annals of Discrete Mathematics (33) Proceedings of the International Conference on Finite 
Geometries and Combinatorial Structures, Vol. 144 of North-Holland Mathematics Studies, North-Holland, 1987, pp. 
41-57. doi:10.1016/30304-0208(08)73047-0 

[31] M. Caselle, F. Gliozzi, S. Necco, Thermal operators and cluster topology in the q -state potts model, J. Phys. A 34 (3) 
(2001) 351. 

[32] G. Grimmett, The stochastic Random-Cluster process and the uniqueness of Random-Cluster measures, Ann. Probab. 
23 (4). 

[33] C.-K. Hu, J.-A. Chen, N. S. Izmailian, P. Kleban, Geometry, thermodynamics, and finite-size corrections in the critical 
Potts model, Phys. Rev. E 60 (6) (1999) 6491. 

[34] A. E. Ferdinand, M. E. Fisher, Bounded and inhomogeneous Ising models. 1. Specific heat anomaly of a finite lattice, 
Phys. Rev. 185 (1969) 832. doi:10.1103/PhysRev.185.832 

[35] J. Salas, A. D. Sokal, Logarithmic corrections and finite-size scaling in the two-dimensional 4-state potts model, J. 
Stat. Phys. 88 (1997) 567. 


29 


[36] M. Nauenberg, D. J. Scalapino, Singularities and Scaling Functions at the Potts-Model Multicritical Point, Phys. Rev. 
Lett. 44 (1980) 837-840. doi:10.1103/PhysRevLett.44.837 

[37] J. L. Cardy, M. Nauenberg, D. J. Scalapino, Scaling theory of the Potts-model multicritical point, Phys. Rev. B 22 
(1980) 2560-2568. doi:10.1103/PhysRevB.22.2560 

[38] Y. Deng, X. Yang, Finite-size scaling of energylike quantities in percolation, Phys. Rev. E 73 (6) (2006) 066116. 

[39] J. M. Schmidt, A simple test on 2-vertex-and 2-edge-connectivity, Inform. Process. Lett. 113 (7) (2013) 241-244. 

[40] Y. Deng, X.-W. Liu, J. L. Jacobsen, Recursive percolation, arXiv preprint arXiv: 1410.3603 

[41] V. Privman, Finite-Size scaling theory, in: V. Privman (Ed.), Finite Size Scaling and Numerical Simulation of Statistical 
Systems, World Scientific, Singapore, 1990, pp. 1-98. 

[42] J. Salas, A. D. Sokal, Universal amplitude ratios in the critical two-dimensional ising model on a torus, J. Stat. Phys. 
98 (3-4) (2000) 551-588, arXiv:9904038vl. 

[43] W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery, Numerical Recipes: The Art of Scientific Computing, 
3rd Edition, Cambridge University Press, Cambridge, 2007. 

[44] A. Coniglio, Cluster structure near the percolation transition, J. Phys. A 15 (1982) 3829. 

[45] B. Duplantier, H. Saleur, Exact fractal dimension of of 2D ising clusters, Phys. Rev. Lett. 63 (1989) 2536. 

[46] Y. Deng, W. Zhang, T. M. Garoni, A. D. Sokal, A. Sportiello, Some geometric critical exponents for percolation and 
the random-cluster model, Phys. Rev. E 81 (2010) 020102. doi: 10.1103/physreve. 81.020102 

[47] V. Beffara, P. Nolin, On monochromatic arm exponents for 2D critical percolation, Ann. Probab. 39 (2011) 1286-1304. 

[48] R. Vasseur, J. L. Jacobsen, H. Saleur, Logarithmic observables in critical percolation, J. Stat. Mech.: Theory Exp. 
2012 (07) (2012) L07001. doi:10.1088/1742-5468/2012/07/L07001 

[49] R. Vasseur, J. L. Jacobsen, Operator content of the critical Potts model in d dimensions and logarithmic correlations, 
Nucl. Phys. B 880 (2014) 435. 

[50] P. Grassberger, Conductivity exponent and backbone dimension in 2-d percolation, Physica A 262 (3) (1999) 251—263. 
doi:10.1016/S0378-4371(98)00435-X 

[51] E. Domany, D. Mukamel, B. Nienhuis, A. Schwimmer, Duality relations and equivalences for models with 0(N) and 
cubic symmetry, Nucl. Phys. B 190 (1981) 279. 

[52] Y. Deng, T. M. Garoni, W. Guo, H. W. J. Blote, A. D. Sokal, Cluster Simulations of Loop Models on Two-Dimensional 
Lattices, Phys. Rev. Lett. 98 (12) (2007) 120601. doi: 10.1103/PhysRevLett. 98.120601 


30 


